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Abstract 

An analytical method of electromagnetic wave interactions with a general radially anisotropic 
cloak is established. It is able to deal with arbitrary parameters {srir), fJ^ir), Stir) and fj-t{r)) of a 
radially anisotropic inhomogeneous shell. The general cloaking condition is proposed from the wave 
relations for the first time. We derive the parameters of a novel class of spherical nonlinear cloaks 
and examine its invisibility performance by the proposed method at various nonlinear situations. 
Spherical metamaterial cloaks with improved invisibility performance is achieved with optimal 
nonlinearity in transformation and core-shell ratio. 

PACS numbers: 41.20.Jb, 42.25.Gy, 42.79.Dj 
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I. INTRODUCTION 



Coordinate transformation for the design process of the cloaking devices has 

received great attention. The cyhndrical/spherical cloaking idea proposed by Pendry [1] is to 
employ radial anisotropic materials whose parameters are determined from the topological 
variation between the original and transformed spaces, based on the invariance of Maxwell's 
equations throughout a specific coordinate transformation J3|]. The idea of cylindrical cloak- 
ing was confirmed by analytical/full- wave methods [s, 6, J and verified by an experiment 



using artificial metamaterials with inclusions of metallic split-ring resonators (SRRs) 
So far, significant progress has been made on the study of cylindrical invisibility cloaks. It 
reveals that the simplified parameters for cylindrical cloaking still allow wave interactions 
with the cloaked object [sl] and the invisibility performance of a cylindrical cloak is very 
sensitive to the geometrical perturbation of its interior boundary 10|], w. 
fixed by introducing PEC/PMC linings onto the inner surface of the shell 



lich can be both 
4, 



Since it is 

challenging to synthesize the magnetic response in optical regime, nonmagnetic cylindrical 
cloaks have been proposed by using quadratic transformation 1^ and the general high-order 
transformation for nonmagnetic cylindrical cloaks in optical frequency is addressed more re- 
cently Nevertheless, it is still difficult to realize the position-dependent cylindrical 
cloak due to the limited resource of natural materials exhibiting radial anisotropy [14^]. In 
view of this, Cai et al. proposed a multilayered cylindrical cloak by dividing the original 
position-dependent cloak into many thin coatings in which the material parameters become 
homogeneous Furthermore, the cylindrical cloak has been theoretically realized by a 
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Arbitrary-shaped 2D cloaks 



19|. 



concentric cylinder of isotropic homogeneous multilayers 
have been investigated theoretically and numerically 171 . 

However, for spherical invisibility cloaks there are still a lot of unknowns to be explored 
because of the complexity in analysis and simulation of scattering properties. Anisotropic 
and position-dependent ideal spherical cloak based on the linear transformation was sug- 
gested by Pendry [l], and it has been shown that spherical cloaks are less sensitive to the 
perturbation than cylindrical cloaks [lo|, which is mathematically proved [20]. There are 
several main streams of studying linear first-order spherical invisibility cloaks, whose the 
required materials and the corresponding methods are distinct. The first approach is the 
classic cloak [1], which is linear, anisotropic and inhomogeneous. In this connection, explic- 
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itly electromagnetic fields have been formulated 



2l| and it is further confirmed that the 



22| . The second is to utilize a homogeneous 



wave cannot interact with the concealed object 
anisotropic metamaterial cover to achieve electromagnetic invisibility [23| via the core-shell 
system. The third is the implem entation of isotropic plasmonic materials as the cloak based 
on cancellation scheme 
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25( 1 ■ The fourth is to substitute the Pendry's classic cloak with 



alternating thin multi-shells and each shell is homogeneous and isotropic 26j. Each approach 
mentioned above has its own advantages and restrictions. For instance, the first approach 
|l| requires higher complexity in material parameters, and the analysis is situated towards 
a particular anisotropy ratio, which is addressed in 26[ . The second approach 23|] removes 
the requirement of material inhomogeneity, in which parameters are position-independent. 
However its cloaking property is quite reliant on the core-shell ratio, and the same feature is 
possessed by the third approach (cancellation scheme). The fourth method has less restric- 
tions on materials but needs a lot of coatings which are sufficiently thin compared with the 
wavelength. The high-order term in the refractive index of an inhomogeneous spherical lens 
is discussed and its possibility of realizing a spherical cloak without parametric singularity 
is addressed [27]. The critical material singularity is thus transformed into the geometrical 



singularity which is less demanding 



28|. 



In this paper, a more general high-order nonlinear transformation will be considered for 
spherical cloaks. We first propose a general algorithm to study the electromagnetic scattering 
by a particle coated by a radially inhomogeneous shell whose anisotropic parameters can be 
arbitrary. We discretize the shell into multiple spherical shells, each of which is homogeneous 
and anisotropic. Also, we propose a novel class of nonlinear transformation based spherical 
cloaks, whose anisotropy ratio is position-dependent and also too complicated to be treated 
by any mentioned methods. By utilizing the established general scattering algorithm, the 
invisibility performance and its dependence on the nonlinear transformation are investigated. 
Finally, the numerical results suggest a particular type of nonlinear spherical cloak providing 
better invisibility than Pendry's linear spherical cloak. 
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FIG. 1: The geometry of the spherical cloak structure. Incident plane wave is propagating along 
z-direction and its electric field is polarized along x-direction. The supscripts and 1 denote the 
parameters of the host and cloaked media, respectively. The anisotropic e and /I presents the 
parameters of the cloak shell. 

II. SCATTERING ALGORITHM FOR A GENERAL RADIALLY ANIOSTROPIC 
METAMATERIAL CLOAK 

Fig. [1] illustrates the configuration of the cloak structure, i.e. the inner and outer radii 
are denoted by a and h respectively; innermost region is filled by an isotropic dielectric 
material characterized by e^^^ and /x*-^-*; intermediate region is occupied by a general spherical 
metamaterial cloak characterized by e and Jl 

£ = er{r)er ^ er + et{r)It, Jl = fir{r)er er + fit{r)It, (1) 

where Er and fir are the radial permittivity and permeability, et and fit are the transversal 
material parameters. It = I — Gr ^ Br = eg ® eg + e^p (E) e^^ is the projection operator onto 
the plane perpendicular to the vector e^, / is the unit three-dimensional dyad, unit vectors 
e^, eg, and e^ are the basis vectors of the spherical coordinates. 

In this section, the scattering theory of multilayer anisotropic spherical particles is pro- 
vided and applied to study a cloak. We suppose that the arbitrary field distribution of the 
incident monochromatic wave interacts with the two-layer sphere. 

Using the separation of variables, the solution of Maxwell's equations in spherical coor- 
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dinates (r, 5, ip) can be presented as 

E(r, e, if) = FiUO, ¥^)E(r), H(r, 9, if) = Fi^O, </^)H(r), (2) 

where the designation E(r) means that the components of the electric field vector depend 
only on the radial coordinate r as Er{r), Eg{r), and E^f,{r) (however, the vector itself includes 
the angle dependence in the basis vectors), and the second rank tensor in three-dimensional 
space Fim serves to separate the variables (/ and m are the integer numbers). It can be 
written as the sum of dyads: 



Im 



YimGr ®er + X;^ ® Be + (e^ X X;^) (g) e^. (3) 



where Yim{6,ip) and X;m(6', </)) are the scalar and vector spherical harmonics, the orthogo- 



nality of which has been well described in SOj] . Tensor functions Fim are very useful, because 



they completely describe the angle dependence of the spherical electromagnetic waves and 
satisfy the orthogonality conditions 

Jo Jo 

where the superscript + stands for the Hermitian conjugate. 

From the commutation of e, JI and Fim, it follows that the electric and magnetic fields 
obey the set of ordinary differential equations 

^dH , 1 i\A(^TT) 



dE 1 iV^(/ + l) 
— -I- -e^ E 

dr r r 



< — + I<E - ^- e^E = iA;o7I ■ H, (5) 



where /cq = cu/c is the wavenumber in vacuum, and uj denotes the circular frequency of the 



incident electromagnetic wave. Quantity is called the tensor dual to the vector n [31 1. 
It results in the vector product, if multiplied by a vector a as n^a = na^ = n x a. 

Eq. ([5]) results from variable separation in Maxwell's equations, and includes two algebraic 
scalar equations, therefore, two field components, Hf and E^., can be expressed by means of 
the rest four components. This can be presented in terms of the matrix link between the 
total fields H = Ht + iJ,.e,, and E = Et + E^e^ and their tangential components Ht and Eti 



( H(r) \ 


( Ht(r) \ 


= V{r) 


\ Et(r) ) 


\ E(r) ) 



(6) 



Excluding the radial components of the fields from Eq. ([5]), we arrive at a set of ordinary 
differential equations of the first order for the tangential components, which can be inter- 
connected by a four- dimensional vector W(r) as 

dWfr 



dr 



ifcoM(r)W(r) 



(7) 



where 



M 




A = D 



w 



-e^ (g) eg, C = -^it{r)el 




lil + l) 



Since Ht and Et are continuous at the spherical interface, they can be used for solving 
the scattering problem. Now, we analyze the situation of r-dependent permittivities and 
permeabilities, which arise from the spherical cloaking. Excluding the (/^-components of 
fields from Eq. ([7j), we derive the differential equation of the second order for the vector 
We = ■ W = {He,Ee): 





St 

^ 



et 

^ 



et. 

£r 



^ 



We = 0, 



(9) 



where the prime denotes the derivative with respect to r. Further we will apply one condition 
on the medium parameters, which is usually used for the spherical cloaks: er{r) = ^r{r) and 
St{r) = due to the impedance matching. Then the equations for and Eg coincide 

and can be written in the form 



,22 _ ^ _ 



e[ l{l + l)et 



We = 0. 



(10) 



ret £r 

This equation can be solved analytically only in very few cases. As an example, we can 
offer a case of Et = fit = cii/^ and Sr = fir = 0'2/^-, where ai^2 are arbitrary values. However 
such dependencies do not provide the cloak properties. Another solvable case in Eq. ( fTOl) is 
just Pendry's cloak, that is, et = fit = b/{h — a) and Sr = fJ^^tir — a)^/r^. 

Although analytical solutions cannot be found for all situations, the general structure 
of solutions can be studied. The solution of two differential equations of the second order 



Eq. ([9]) contains four integration constants Ci, C2, c'l, and C2. The constants can be joined 
together into a couple of vectors Ci = ciBq + c[e^ and C2 = C2eg + Cge^. (/^-components of the 
field vectors and are expressed in terms of the 6'-components which have been already 
determined. The relation between 6- and y9-components follows from Eq. ([7]). Summing up 
both components, the resultant field can be presented as 

W = S{r)C, S{r) = I "^'^'^ "^'^'^ I , C = I , (11) 




Ci(r) C2(r) 

where 7712 and Ci,2 are the two-dimensional blocks of the matrix S{r). {rji, (i, Ci} and 
{V2, C2, C2} denote the first and second sets of the independent solution of Eq. ([9]), re- 
spectively. Therefore, the general solution can be decomposed into two terms as W = 
WW + W(2), where 

W« = T" I = I "■ I c. W'^' = = I * I c. (12) 







Electric and magnetic fields of each independent wave are related by means of impedance 
tensor T as Etj = Tjiitj (j = 1,2). Thus, the impedance tensor equals 

T,ir) = Qir)r^T\r). (13) 

Vectors Ci and C2 can be expressed by means of the known tangential electromagnetic 
field W(a) as C = S~^{a)W{a). Thus Eq. ffTTl) can be rewritten as follows 

W(r)=fi^W(a), n: = S{r)S-\a), (14) 

where the evolution operator (transfer matrix) fi^ connects tangential field components 
at two distinct spatial points, i.e., r and a. One can obtain the complete solution of the 
fields E(r) and H(r) by summing over I and m in the subsequent tensor product of Fim{6, ^p) 
describing angle dependence (Eq. ([3])), matrix y'(r) restoring the fields with their tangential 
components (Eq. ([6])), and tangential field vectors (Eq. ffTTl) ): 




H(r) \ ^ ^ FUd.v) \ / 7^\{r) r/^(r) 

V [r) 



\ ^2 



(15) 



In general, the solutions cannot be studied in the closed form for nonlinear spherical 
cloaks. Therefore, the approximate method of numerical computations is applied. An 



inhomogeneous anisotropic spherical cloak a < r < b is equally divided into homogeneous 
anisotropic spherical layers, i.e., replaced with a multi-layer structure. The number of the 
layers strongly determines the accuracy of calculations. The j-th homogeneous shell is 
situated in the region between r = aj_i and r = aj, where j = 1,...,N, = a and 
= b. Wave solution of the single homogeneous layer can be represented in the form of 
the evolution operator Qa^^i- The solution for the whole inhomogeneous shell is thus the 
subsequent product of the elementary evolution operators: 

nt = nl ...fi^^fi^^ (16) 

The solution of Eq. (Q in one layer with constant permittivities Er, St and permeabilities 
/ir, /it is expressed by means of a couple of independent spherical functions g^"* and g^^: 

9i^Xktr)ci + gif{ktr)c2 
gl^JihrK + gifiktry^ 





H 







(17) 



where /c+ 



ko,/e^t, 1^1 = yi(]TTK7^7Tl74 - 1/2, z/2 = 



Wl^+ 1/4 - 1/2 



which applies to both uniaxial anisotropic and bianisotropic media [281 . 1321 . l33| . Functions 
gl^'"^^ of the order u can be spherical Bessel functions, modified spherical Bessel functions, 
or spherical Hankel functions. 

Blocks 1] and ( introduced in Eq. ffTTl) are the tensors 



111,2 
Cl,2 



.(1,2) 



eg (g) eg 



d{rg^ 



(1,2)^ 
1/2 ) 



Htkor 
i 



d(r^^| 



dr ^ 



eg, 



(18) 



Efkor dr 

which can be presented as two-dimensional matrices for computation purposes. 

Now we turn to the scattering of electromagnetic waves from the cloaking structure 
depicted in Fig. [1] We suppose that an electromagnetic field Hinc(r) and Einc(r) is incident 
on the coated spherical particle from air {e^^^ = 1, = 1). Wave solutions in each of the N 
layers can be written using the general solution Eq. f[T^ . which is already known. Scattered 
field propagating in air can be presented by the superposition of diverging spherical waves 
which are mathematically described by spherical Hankel functions of the first kind hl^\x). 
Let us first introduce rj and ( which correspond to the tensors rj and ( in Eq. (1181) when 
Hankel functions replace Bessel functions. Then we obtain the scattered fields 



Es, 



[r) 




rHr) 



fj\r){fj\b))-'Hi^{b), 



(19) 



where = C'(^0 ^ is the impedance tensor of the Ith scattered wave, and H^^(6) is the 
tangential magnetic field at the outer interface r = b. Applying the evolution operator fi^, 
the electromagnetic field in the shell takes the form 









( 



Im 



Hf(a), 



[T[{a) 



(20) 



where T[ = CU^i) ^ impedance tensor of the Ith wave inside the inner sphere (region 



1), and H]™'(a) is the tangential magnetic field at the inner interface of the shell r 



a. 



By projecting the fields onto the outer interface r = b and integrating over the angles 9 
and if with the help of orthogonality condition Eq. (jl]) , we derive the boundary conditions 



W 



Im 



r'(6) ) I T[{a) 



(21) 



where 



W 



Im 



Jo 



sin 9d9dip. 



(22) 



Eq. fl2ip represents the system of four linear equations for four components of the vectors 
Hg™ and H']™. Finally, one can derive the amplitude of the scattered electromagnetic field 
(see Appendix) 

-1 



ui^ib) 



r[ia) -I Hit 



r\b) 



T[ia) -I 



jlm 
inc 



(23) 



The scattered field in far zone can be characterized by the differential cross-section (power 
radiated to e^-direction per solid angle do) 



da 



|Hse(r) 



(24) 



do |Hinc(r)r 

In our notations, the differential cross-section averaged over the azimuthal angle if (over 
polarizations) normalized by the geometrical cross-section (jg = vrfe^ takes the form 

2 



da 



1 



an sin 9d9 cr„|Hinr-P 



E 



mc| m=-oo 



Y: r'-'FU0,o)r,r\b)iii^{b) 



l=\m\ 



(25) 



From the point of view of the scattering theory, it is straightforward to define a cloak 
as one specially matched layer that provides zero scattering for arbitrary materials inside. 
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In [21j zero scattering was proved analytically for the Pendry's spherical cloak. Here, we 
have proposed a more general scattering algorithm for radially anisotropic materials, which 
is useful in studying the scattering of spherical cloaks based on complex (e.g., high-order, 
nonlinear, etc.) transformations. From the proposed scattering theorem, we can determine 
the invisibility condition (zero scattering) specified by the condition Hg™(6) = 0, which in 
turn can be rewritten using Eq. (!23|) as follows 

r[{a) -I^Qty^t = 0. (26) 

Arbitrary incident electromagnetic field W(j^ can be excluded from this expression. In fact, 
zero scattering can be obtained for the trivial situation: electromagnetic field is scattered by 
a "virtual" air sphere at radius b. This assumption can be presented in the form analogous 
to Eq. (EE]): 

V|,(6) -/)w|- = 0, (27) 

where Fq is the impedance tensor of the Ith wave in the "virtual" air sphere. Hence, we 
have the relation 



r{{a) -l)^t=[P,{b) -I). (28) 

Impedance tensor F^^ contains material parameters e^^^ and /^'-^^ of the sphere inside the 
cloaking shell. At the same time, zero scattering should be held for arbitrary e^^^ and 
jj,^^^ of the inner core. It implies that the partial derivative of Eq. fl28|) with respect to 
e^^^ needs to be zero for arbitrary e^^^ to satisfy the zero scattering condition. Note that 
only the impedance tensor F'^^ contains e^^\ and therefore the right-hand side of Eq. fl25]) 

vanishes after the differentiation, which results in — ^ — —Q^; + ( r[{a) —I ) = 

l^i- ) fife = -^-(T) [ I ) fifo = 0. It is now straightforward that we need the relation to 



be satisfied, i.e., ( / ) = 0. 



By multiplying this equation by T[{a) and subtracting it from Eq. fl28|) . we arrive at the 
equation: ^0 / ^ = ^ — Fq(6) / ^, which does not contain the material parameters of 
the inner sphere. Finally, we derive the evolution operator of the cloaking layer: 

^t=i : , " I • (29) 




10 



This condition defines the cloak and can be satisfied for the specially chosen evolution 
operator of the cloaking shell. The evolution operator is the degenerate block matrix, 
whose inverse matrix is not defined. It should be noted that relation Eq. fl29l) is indepen- 
dent on the material of the inner sphere. On the other hand, the derived relation connects 
the wave solutions in the cloak (evolution operator Q^) with wave solutions in the "equiva- 
lent" homogeneous air sphere (impedance tensor Fq). Therefore, it effectively performs the 
coordinate transformation for the solutions, but not for the material parameters as usual. 

Substitution of the cloaking condition H'™(6) = into Eq. (12T]) results in 



From Eq. ([29D, it becomes clear that li^{^{a) = and E^^{a) = r{{a)H^{^{a) = 0. Thus one 
may conclude that both electric and magnetic fields equal zero at the boundary r = a, and 
therefore there is no field at any spatial point inside the inner sphere. In the cloaking shell, 
the fields equal zero at the inner boundary r = a (owing to the continuity of the tangential 
fields) and equal incident fields at the outer boundary r = b. Then, the field inside the cloak 
must obey (see Eq. (!20|) ) 



III. NONLINEAR TRANSFORMATION BASED SPHERICAL CLOAKS 

Now, let us consider a novel class of the nonlinear transformation based (NTB) spherical 
cloak, whose EM interaction can be characterized by the proposed scattering theorem. Fig. [1] 
can be regarded as the compressed space (r) from the original space (r'), i.e., an air sphere 
< r' < 6. We propose a nonlinear transformation function 



which obviously satisfies the transformation (when r = a, r' = 0; and when r = b, r' = b). 
The value of "x" is a factor to control the nonlinearity degree in the transformation, which 
can be arbitrary from to oo. 

Due to the invariance of Maxwell's equation under coordinate transformations from the 
original space to transformed space, the parameters (e, Jl) in the shell of Fig. [1] can be 




(30) 




(31) 




(32) 
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expressed in terms of those parameters in the original space, i.e., s' = 1 and ^t' = 1, 

£ = AA^I det ( A) , 71 = AA^I det {A) , (33) 

where A is the Jacobian matrix with elements Aij = dri/dr'-. 

One can see that the proposed prescribed function Eq. fl32|) is only dependent on radial 
position r. Then it is easy to find that the Jacobian matrix is diagonal, and Eq. (|33l) can 
thus be rewritten as 

£ = 71 = diag[A^, A^, Ajl/Ar-AeA^ = diag[-^, -^], (34) 

M'^cj) '^r^(f> \M 

where the principal stretches of the Jacobian matrix are 

^ dr {b — aYr^^^ ^ ^ r [b — aYr^^^ 
^ dr' xa6^+^(r — a)^^^ ' ^ j-i i)X+i(j, _ 

Finally, one can obtain the parameters of the NTB cloak (a < r < 6) in Fig. [T] 

6^+^(r - a)^+^ 



fir 



xa{b — a)^r 



eg= fi9 = e^ = _ ^y^^^+i ■ (36) 

Such ideal NTB spherical cloak is difficult to be fabricated in practice. However, to some 
extent, it can be alleviated by dividing the inhomogeneous cloak shell into homogeneous 
multilayers. The case for cylindrical cloaks has been studied and it shows that only several 



optimized layers can achieve the invisibility 



34( 1 . Here, the optimization is out of the scope 



of this paper. Our paper is to reveal some novel NTB spherical cloaks which provides 
better invisibility performance than Pendry's classic one, based on the proposed general 
scattering theory. The realistic NTB cloaks can be produced using sputtering techniques, so 
that a number of discrete layers should be applied over the spherical core. To demonstrate 
the capability of the proposed spherical cloaks, we present the differential cross-sections 
normalized by the geometric cross-section of the cloak (see Eq. (1251) ). 

In Fig. [21 we analyze the dependence on the total number of the layers dividing the 
cloaking shell. The increase in the layer number gives rise to more accurate approximation of 
the original inhomogeneous model in Fig. [1], and the decrease of the scattering cross-section 
is expected with the increase of the number A^. If one uses the present scattering method 
with = 50 to divide both Pendry's linear cloak and a specific NTB cloak at x = 1, 
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the forward scattering is approximately the same, but over the whole range of scattering 
angles Pendry's cloak presents better invisibility. The following discussion will address the 
importance of this nonlinear factor x in beating the classic linear cloak. 




Scattering Angie e (degree) 

FIG. 2: Differential cross-section da /{ag sin 9 d9) of the NTB spherical cloak (x = 1) for different 
number of the layers dividing the inhomogeneous coating. Parameters: e^^^ = 1.45^, fj,^-^^ = 1.0, 
koa = vr, kQb = 2tt. 

The parameter x is a convenient tool to control the quality of the NTB spherical cloak. 
We assume the number of homogeneous sublayers = 50 for all following simulations. 
In Fig. [3] the differential cross-sections at different x are demonstrated. If x is less than 
unity, the cross-section is inversely proportional to x, which is not desired in the sense of 
invisibility. For NTB spherical cloaks with x < 1, the cloaking performance is degraded due 
to the abrupt increase of the transverse dielectric permittivity et near the inner interface 
r = a of the clad (see Fig. HI). 

The abrupt change of the material parameters is undesirable not only for invisibility 
performance but also for the realization point of view. The radial dielectric permittivity 
behaves in a similarly monotonic way for all values of x. The dependence appears to be 
mainly linear for small parameters x. The dependence of the transverse dielectric permittiv- 
ity is more complicated. One particular NTB cloak is realized at the parameter x = 2 when 
the transverse permittivity in the cloak becomes non-monotonic and eventually returns to 
Et = 2 aX r = b, which provides even lower cross-section over whole observation angles than 
the Pendry's cloak does. If we compare the dielectric permittivities of the proposed NTB 
cloak with that of Pendry's cloak (Fig. H]), it can be noted that the dependence of radial 
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Scattering Angle e (degree) 

FIG. 3: Differential cross-section of tlie nonlinear cloak with different parameters x. Parameters: 
e(i) = 1.452, = 1.0, koa = vr, kob = 2-k,N = 50. 




Clad Thickness k d 



FIG. 4: Radial Sr and transverse ej dielectric permittivities of the cloaking shell with different 
parameters x. Cloak is extended from koa = vr to kob = 2tt. 

permittivity are still quite close to each other in the cloak region. However, one may 
ask whether a; = 2 is the only choice or not. In Fig. [5], it gives the answer that in the 
sense of total cross section, there is a range of x in which the proposed NTB spherical cloak 
outperforms the classic linear spherical cloak. When "x" increases and jumps out of this 
optimal region, the cloaking effects compared with Pendry's cloak are degraded which can 
be verified by the bistatic cross-section of x = 1 and a; = 3 in Fig. [31 

Furthermore, we investigate how the transverse permittivity St varies near the optimal 
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-50-1 




-60-1 . , . , . , . ,— 

1,0 1,5 2,0 2,5 3,0 

Parameter x 

FIG. 5: Scattering cross-section versus x for NTB spherical cloaks. Parameters: e*-"^^ = 1.45^, 
/i^-*^) = 1.0, koa = vr, kob = 2tt, N = 50. The range of x where SCS is lower than that of Pendry's 
spherical cloak is the optimal region of x for desired NTB spherical cloaks. 




FIG. 6: Variance of transverse permittivity £t along the radial direction in the region of the cloaking 
shell {a < r < b) under different values of x near the optimal range as shown in Fig. [5j Parameters: 
k^a = vr and kob = 27r. 

region of "x" in Fig. [51 From Fig. [6|, it can be observed that: 1) when x is slightly above 1, 
the requirement of 6t near the inner boundary r = a drops significantly compared to those 
curves whose "x" is smaller than 1 in Fig. \5\ 2) when "x" becomes larger and larger, the et 
at the outer boundary r = b turns to be more deviated from that value of its corresponding 
Er', 3) when x falls into the optimal region, those curves of transverse permittivities are 
non-monotonic along the radial direction in the cloaking shell, and their maxima and overall 
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values of et are smaller than those whose x becomes further smaller or larger. These explain 
why there exist an optimal region for "x" where the total scattering cross-section can be 
lower than Pendry's classic one. Also, it provides us another way to predict whether a 
certain "x" for a NTB spherical cloak is optimal or not. 

Now we continue to study the dependence of its invisibility upon the ratio b/a of the 
particular NTB spherical cloak with x = 2 which is discretized into = 50 layers. We keep 
the inner radius a unchanged. In Fig. El^a), different ratios of b/a are considered. Compared 
with the other three values of b/a, it seems that b/a = 2 provides the best cloaking effects at 
nearly all angles (except for the angle at 52°) for the x = 2 NTB cloak. Another interesting 
finding is that: when b/a > 2, the cross-section will be higher than that of 6/a = 2 over the 
whole range of angles; when 6/a — > 1, though the angle-averaged cross-section will still be 
higher but at certain angles, its cross-section could be lower than that of b/a = 2. From the 
view of total scattering, it is important to consider how the ratio b/a should be selected for 
X = 2 NTB cloak so as to provide improved cloaking. In Fig. [7](b), one can clearly see the 
optimal domain of b/a in which the cross-section is smaller than Pendry's spherical cloak. 
Certainly, all values x within the desired region for the purpose of improved cloaking in 
Fig. [5] have their corresponding domain of optimal b/a. 




20 40 60 80 100 120 140 160 180 1,5 2,0 2,5 

Scattering Angle e (degree) p^gd^ ^/^ 



FIG. 7: The role of core-shell ratio b/a in the cloaking improvement for x = 2 NTB spherical cloak: 
(a) differential cross-section versus angle at selected ratios; (b) scattering cross-section versus ratio 
b/a. Parameters: e^^^ = 1.45^, ^^^^ = 1.0, koa = vr, and N = 50. 

Let us compare the scattering diagrams for both a bare glass sphere and a cloaked glass 
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sphere. We take the clad in the form of the nonhnear cloak with parameter x = 2. From 
Fig. [HI we see that the cloaking shell noticeably reduces the scattering. Another general 
property of the cloaks, i.e, the exactly diminished backscattering, is also present in the 
figure. In what follows, we consider their respective near-field wave interactions, which 
correspond to far-field results in Fig. [HI 



FIG. 8: Differential cross-section of a glass sphere with and without cloaking shell. Parameters: 
e(i) = 1.452, ^(1) = 1.0, koa = vr, kob = 2tt, x = 2, N = 50. 

The near-field perturbation of the cloaked and non-cloaked glass particles is demonstrated 
in Fig.[9l Comparing Fig.[9]^a) with Fig.[9](b), the invisibility performance is well pronounced. 
In Fig. [9](a) the EM wave travels only through the clad and takes near zero values in the 
vicinity of the inner radius a. The field does not enter the glass core, therefore it does not 
matter that which material is situated therein. If the cloak is less ideal than that shown in 
the figure, the incident field will be scattered by the spherical particle and will penetrate 
the glass core, i.e., the object becomes visible. 

IV. CONCLUSION 

We have studied the manifold of the nonlinear cloaks differing in parameter x. Since 
there is no closed form solution for the proposed nonlinear cloaks, an approximate model by 
replacing an inhomogeneous shell with homogeneous spherical layers has been numerically 
analyzed, with the help of the proposed scattering algorithm for multilayered rotationally 
anisotropic shells. The general cloaking condition was derived from the scattering algorithm, 
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FIG. 9: Real part of the electric field on x-z plane scattered (a) by the cloaking shell gathered 
round the glass core and (b) by the glass core itself. Parameters: e^^^ = 1.45^, n^^^ = 1.0, koa = tt, 
kob = 27r, X = 2, iV = 50. 

which is in contrast to the method of coordinate transform. We have also demonstrated that 
better approximate spherical cloaks can be realized by properly choosing parameters x and 
b/a. In practical applications, such a class of NTB spherical cloaks can provide improved 
invisibility performance. 

Appendix A 



In order to exclude the constant vector Hi , one should multiply Eq. (EI]) by fi^ = {^1) 



b\-l. 



r'(6) ) I T[ia) 



Hf (a) 



(A-1) 



Then Eq. flA-ll) is further multiplied by the block matrix | T[{a) —I ), and the right-hand 
side vanishes: 

/ 



T[{a) -/)^^^W|-+(rUa) 



iiim = 0. 



(A-2) 



On the other hand, vector H^™ can be obtained in a similar way. 
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Abstract 

An analytical method of electromagnetic wave interactions with a general radially anisotropic 
cloak is established. It is able to deal with arbitrary parameters {srir), fJ^ir), Stir) and fj-t{r)) of a 
radially anisotropic inhomogeneous shell. The general cloaking condition is proposed from the wave 
relations for the first time. We derive the parameters of a novel class of spherical nonlinear cloaks 
and examine its invisibility performance by the proposed method at various nonlinear situations. 
Spherical metamaterial cloaks with improved invisibility performance is achieved with optimal 
nonlinearity in transformation and core-shell ratio. 

PACS numbers: 41.20.Jb, 42.25.Gy, 42.79.Dj 



* Electronic address: 



cwq@mit.edu 



I. INTRODUCTION 



mm 



Coordinate transformation for the design process of the cloaking devices has 

received great attention. The cyhndrical/spherical cloaking idea proposed by Pendry [1] is to 
employ radial anisotropic materials whose parameters are determined from the topological 
variation between the original and transformed spaces, based on the invariance of Maxwell's 
equations throughout a specific coordinate transformation J3|]. The idea of cylindrical cloak- 
ing was confirmed by analytical/full- wave methods [s, 6, J and verified by an experiment 



using artificial metamaterials with inclusions of metallic split-ring resonators (SRRs) 
So far, significant progress has been made on the study of cylindrical invisibility cloaks. It 
reveals that the simplified parameters for cylindrical cloaking still allow wave interactions 
with the cloaked object 9i] and the invisibility performance of a cylindrical cloak is very 



sensitive to the geometrical perturbation of its interior boundary 10|], w. 
fixed by introducing PEC/PMC linings onto the inner surface of the shell 



lich can be both 
4, 



Since it is 

challenging to synthesize the magnetic response in optical regime, nonmagnetic cylindrical 
cloaks have been proposed by using quadratic transformation 12| and the general high- 



order transformation for nonmagnetic cylindrical cloaks in optical frequency is addressed 
more recently [isl. Nevertheless, it is still difficult to realize the position-depedent cylindri- 
cal cloak due to the limited resource of natural materials exhibiting radial anisotropy |l6 |. 
In view of this, Cai et al. proposed a multilayered cylindrical cloak by dividing the original 
position-dependent cloak into many thin coatings in which the material parameters become 
homogeneous 17|]. Furthermore, the cylindrical cloak has been theoretically realized by a 
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2d, 



Arbitrary-shaped 2D cloaks 



2l|. 



concentric cylinder of isotropic homogeneous multilayers 
have been investigated theoretically and numerically 19|, 

However, for spherical invisibility cloaks there are still a lot of unknowns to be explored 
because of the complexity in analysis and simulation of scattering properties. Anisotropic 
and position-dependent ideal spherical cloak based on the linear transformation was sug- 
gested by Pendry [l], and it has been shown that spherical cloaks are less sensitive to the 
perturbation than cylindrical cloaks [l^ . There are several main streams of studying linear 
first-order spherical invisibility cloaks, whose the required materials and the corresponding 
methods are distinct. The first approach is the classic cloak [1], which is linear, anisotropic 
and inhomogeneous. In this connection, explicitly electromagnetic fields have been formu- 



2 



ated 



22| and it is further confirmed that the wave cannot interact with the concealed object 



231] . The second is to utihze a homogeneous anisotropic metamaterial cover to achieve elec- 



tromagnetic invisibility 26| via the core-shell system. The third is the implementation of 



isotropic plasmonic materials as the cloak based on cancellation scheme |24j. |25[| . The fourth 
is to substitute the Pendry's classic cloak with alternating thin multi-shells and each shell 
is homogeneous and isotropic 27|. Each approach mentioned above has its own advantages 
and restrictions. For instance, the first approach |l| requires higher complexity in mate- 
rial parameters, and the analysis is situated towards a particular anisotropy ratio, which 
is addressed in The second approach 26] removes the requirement of material inho- 
mogeneity, which parameters are position-independent. However its cloaking property is 
quite reliant on the core-shell ratio, and the same feature is possessed by the third approach 
(cancellation scheme). The fourth method has less restrictions on materials but needs a lot 
of coatings which are sufficiently thin compared with the wavelength. The high-order term 
in the refractive index of an inhomogeneous spherical lens is discussed and its possibility of 
realizing a spherical cloak without parametric singularity is addressed [l^ . The critical ma- 
;erial singularity is thus tranformed into the geometrical singularity which is less demanding 

In this paper, a more general high-order nonlinear transformation will be considered for 
spherical cloaks. We first propose a general algorithm to study the electromagnetic scattering 
by a particle coated by a radially inhomogeneous shell whose anisotropic parameters can be 
arbitrary. We discretize the shell into multiple spherical shells, each of which is homogeneous 
and anisotropic. Also, we propose a novel class of nonlinear transformation based spherical 
cloaks, whose anisotropy ratio is position-dependent and also too complicated to be treated 
by any mentioned methods. By utilizing the established general scattering algorithm, the 
invisibility performance and its dependence on the nonlinear transformation are investigated. 
Finally, the numerical results suggest a particular type of nonlinear spherical cloak providing 
better invisibility than Pendry's linear spherical cloak. 
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FIG. 1: The geometry of the spherical cloak structure. Incident plane wave is propagating along 
z-direction and its electric field is polarized along x-direction. The supscripts and 1 denote the 
parameters of the host and cloaked media, respectively. The anisotropic e and /I presents the 
parameters of the cloak shell. 

II. SCATTERING ALGORITHM FOR A GENERAL RADIALLY ANIOSTROPIC 
METAMATERIAL CLOAK 

Fig. [1] illustrates the configuration of the cloak structure, i.e. the inner and outer radii 
are denoted by a and h respectively; innermost region is filled by an isotropic dielectric 
material characterized by e^^^ and /x*-^-*; intermediate region is occupied by a general spherical 
metamaterial cloak characterized by e and Jl 

£ = er{r)er ^ er + et{r)It, Jl = fir{r)er er + fit{r)It, (1) 

where Er and fir are the radial permittivity and permeability, et and fit are the transversal 
material parameters. It = I — Gr ^ Br = eg ® eg + e^p (E) e^^ is the projection operator onto 
the plane perpendicular to the vector e^, / is the unit three-dimensional dyad, unit vectors 
e^, eg, and e^ are the basis vectors of the spherical coordinates. 

In this section, the scattering theory of multilayer anisotropic spherical particles is pro- 
vided and applied to study a cloak. We suppose that the arbitrary field distribution of the 
incident monochromatic wave interacts with the two-layer sphere. 

Using the separation of variables, the solution of Maxwell's equations in spherical coor- 
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dinates (r, 5, ip) can be presented as 

E(r, e, if) = FiUO, ¥^)E(r), H(r, 9, if) = Fi^O, </^)H(r), (2) 

where the designation E(r) means that the components of the electric field vector depend 
only on the radial coordinate r as Er{r), Eg{r), and E^f,{r) (however, the vector itself includes 
the angle dependence in the basis vectors), and the second rank tensor in three-dimensional 
space Fim serves to separate the variables (/ and m are the integer numbers). It can be 
written as the sum of dyads: 



Im 



YimGr ®er + X;^ ® Be + (e^ X X;^) (g) e^. (3) 



where Yim{6,ip) and X;m(6', </)) are the scalar and vector spherical harmonics, the orthogo- 



nality of which has been well described in 29|] . Tensor functions Fim are very useful, because 



they completely describe the angle dependence of the spherical electromagnetic waves and 
satisfy the orthogonality conditions 

Jo Jo 

where the superscript + stands for the Hermitian conjugate. 

From the commutation of e, JI and Fim, it follows that the electric and magnetic fields 
obey the set of ordinary differential equations 

^dH , 1 i\A(^TT) 



dE 1 iV^(/ + l) 
— -I- -e^ E 

dr r r 



< — + I<E - ^- e^E = iA;o7I ■ H, (5) 



where /cq = cu/c is the wavenumber in vacuum, and uj denotes the circular frequency of the 



incident electromagnetic wave. Quantity is called the tensor dual to the vector n [30|. 
It results in the vector product, if multiplied by a vector a as n^a = na^ = n x a. 

Eq. ([5]) results from variable separation in Maxwell's equations, and includes two algebraic 
scalar equations, therefore, two field components, Hf and E^., can be expressed by means of 
the rest four components. This can be presented in terms of the matrix link between the 
total fields H = Ht + iJ,.e,, and E = Et + E^e^ and their tangential components Ht and Eti 



( H(r) \ 


( Ht(r) \ 


= V{r) 


\ Et(r) ) 


\ E(r) ) 



(6) 



Excluding the radial components of the fields from Eq. ([5]), we arrive at a set of ordinary 
differential equations of the first order for the tangential components, which can be inter- 
connected by a four- dimensional vector W(r) as 

dWfr 



dr 



ifcoM(r)W(r) 



(7) 



where 



M 




A = D 



w 



-e^ (g) eg, C = -^it{r)el 




lil + l) 



Since Ht and Et are continuous at the spherical interface, they can be used for solving 
the scattering problem. Now, we analyze the situation of r-dependent permittivities and 
permeabilities, which arise from the spherical cloaking. Excluding the (/^-components of 
fields from Eq. ([7j), we derive the differential equation of the second order for the vector 
We = ■ W = {He,Ee): 





St 

^ 



et 

^ 



et. 

£r 



^ 



We = 0, 



(9) 



where the prime denotes the derivative with respect to r. Further we will apply one condition 
on the medium parameters, which is usually used for the spherical cloaks: er{r) = ^r{r) and 
St{r) = due to the impedance matching. Then the equations for and Eg coincide 

and can be written in the form 



,22 _ ^ _ 



e[ l{l + l)et 



We = 0. 



(10) 



ret £r 

This equation can be solved analytically only in very few cases. As an example, we can 
offer a case of Et = fit = cii/^ and Sr = fir = 0'2/^-, where ai^2 are arbitrary values. However 
such dependencies do not provide the cloak properties. Another solvable case in Eq. ( fTOl) is 
just Pendry's cloak, that is, et = fit = b/{h — a) and Sr = fJ^^tir — a)^/r^. 

Although analytical solutions cannot be found for all situations, the general structure 
of solutions can be studied. The solution of two differential equations of the second order 



Eq. ([9]) contains four integration constants Ci, C2, c'l, and C2. The constants can be joined 
together into a couple of vectors Ci = ciBq + c[e^ and C2 = C2eg + Cge^. (/^-components of the 
field vectors and are expressed in terms of the 6'-components which have been already 
determined. The relation between 6- and y9-components follows from Eq. ([7]). Summing up 
both components, the resultant field can be presented as 

W = S{r)C, S{r) = I "^'^'^ "^'^'^ I , C = I , (11) 




Ci(r) C2(r) 

where 7712 and Ci,2 are the two-dimensional blocks of the matrix S{r). {rji, (i, Ci} and 
{V2, C2, C2} denote the first and second sets of the independent solution of Eq. ([9]), re- 
spectively. Therefore, the general solution can be decomposed into two terms as W = 
WW + W(2), where 

W« = T" I = I "■ I c. W'^' = = I * I c. (12) 







Electric and magnetic fields of each independent wave are related by means of impedance 
tensor T as Etj = Tjiitj (j = 1,2). Thus, the impedance tensor equals 

T,ir) = Qir)r^T\r). (13) 

Vectors Ci and C2 can be expressed by means of the known tangential electromagnetic 
field W(a) as C = S~^{a)W{a). Thus Eq. ffTTl) can be rewritten as follows 

W(r)=fi^W(a), n: = S{r)S-\a), (14) 

where the evolution operator (transfer matrix) fi^ connects tangential field components 
at two distinct spatial points, i.e., r and a. One can obtain the complete solution of the 
fields E(r) and H(r) by summing over I and m in the subsequent tensor product of Fim{6, ^p) 
describing angle dependence (Eq. ([3])), matrix y'(r) restoring the fields with their tangential 
components (Eq. ([6])), and tangential field vectors (Eq. ffTTl) ): 




H(r) \ ^ ^ FUd.v) \ / 7^\{r) r/^(r) 

V [r) 



\ ^2 



(15) 



In general, the solutions cannot be studied in the closed form for nonlinear spherical 
cloaks. Therefore, the approximate method of numerical computations is applied. An 



inhomogeneous anisotropic spherical cloak a < r < b is equally divided into homogeneous 
anisotropic spherical layers, i.e., replaced with a multi-layer structure. The number of the 
layers strongly determines the accuracy of calculations. The j-th homogeneous shell is 
situated in the region between r = aj_i and r = aj, where j = 1,...,N, = a and 
= b. Wave solution of the single homogeneous layer can be represented in the form of 
the evolution operator Qa^^i- The solution for the whole inhomogeneous shell is thus the 
subsequent product of the elementary evolution operators: 

nt = nl ...fi^^fi^^ (16) 

The solution of Eq. (Q in one layer with constant permittivities Er, St and permeabilities 
/ir, /it is expressed by means of a couple of independent spherical functions g^"* and g^^: 

9i^Xktr)ci + gif{ktr)c2 
gl^JihrK + gifiktry^ 





H 







(17) 



where /c+ 



ko,/e^t, 1^1 = yi(]TTK7^7Tl74 - 1/2, z/2 = 



Wl^+ 1/4 - 1/2 



which applies to both uniaxial anisotropic and bianisotropic media [151 . l3lL l32| . Functions 
gl^'"^^ of the order u can be spherical Bessel functions, modified spherical Bessel functions, 
or spherical Hankel functions. 

Blocks 1] and ( introduced in Eq. ffTTl) are the tensors 



111,2 
Cl,2 



.(1,2) 



eg (g) eg 



d{rg^ 



(1,2)^ 
1/2 ) 



Htkor 
i 



d(r^^| 



dr ^ 



eg, 



(18) 



Efkor dr 

which can be presented as two-dimensional matrices for computation purposes. 

Now we turn to the scattering of electromagnetic waves from the cloaking structure 
depicted in Fig. [1] We suppose that an electromagnetic field Hinc(r) and Einc(r) is incident 
on the coated spherical particle from air {e^^^ = 1, = 1). Wave solutions in each of the N 
layers can be written using the general solution Eq. f[T^ . which is already known. Scattered 
field propagating in air can be presented by the superposition of diverging spherical waves 
which are mathematically described by spherical Hankel functions of the first kind hl^\x). 
Let us first introduce rj and ( which correspond to the tensors rj and ( in Eq. (1181) when 
Hankel functions replace Bessel functions. Then we obtain the scattered fields 



Es, 



[r) 




rHr) 



fj\r){fj\b))-'Hi^{b), 



(19) 



where = C'(^0 ^ is the impedance tensor of the Ith scattered wave, and H^^(6) is the 
tangential magnetic field at the outer interface r = b. Applying the evolution operator fi^, 
the electromagnetic field in the shell takes the form 









( 



Im 



Hf(a), 



[T[{a) 



(20) 



where T[ = CU^i) ^ impedance tensor of the Ith wave inside the inner sphere (region 



1), and H]™'(a) is the tangential magnetic field at the inner interface of the shell r 



a. 



By projecting the fields onto the outer interface r = b and integrating over the angles 9 
and if with the help of orthogonality condition Eq. (jl]) , we derive the boundary conditions 



W 



Im 



r'(6) ) I T[{a) 



(21) 



where 



W 



Im 



Jo 



sin 9d9dip. 



(22) 



Eq. fl2ip represents the system of four linear equations for four components of the vectors 
Hg™ and H']™. Finally, one can derive the amplitude of the scattered electromagnetic field 
(see Appendix) 

-1 



ui^ib) 



r[ia) -I Hit 



r\b) 



T[ia) -I 



jlm 
inc 



(23) 



The scattered field in far zone can be characterized by the differential cross-section (power 
radiated to e^-direction per solid angle do) 



da 



|Hse(r) 



(24) 



do |Hinc(r)r 

In our notations, the differential cross-section averaged over the azimuthal angle if (over 
polarizations) normalized by the geometrical cross-section (jg = vrfe^ takes the form 

2 



da 



1 



an sin 9d9 cr„|Hinr-P 



E 



mc| m=-oo 



Y: r'-'FU0,o)r,r\b)iii^{b) 



l=\m\ 



(25) 



From the point of view of the scattering theory, it is straightforward to define a cloak 
as one specially matched layer that provides zero scattering for arbitrary materials inside. 
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In [22j zero scattering was proved analytically for the Pendry's spherical cloak. Here, we 
have proposed a more general scattering algorithm for radially anisotropic materials, which 
is useful in studying the scattering of spherical cloaks based on complex (e.g., high-order, 
nonlinear, etc.) transformations. From the proposed scattering theorem, we can determine 
the invisibility condition (zero scattering) specified by the condition Hg™(6) = 0, which in 
turn can be rewritten using Eq. (!23|) as follows 

r[{a) -I^Qty^t = 0. (26) 

Arbitrary incident electromagnetic field W(j^ can be excluded from this expression. In fact, 
zero scattering can be obtained for the trivial situation: electromagnetic field is scattered by 
a "virtual" air sphere at radius b. This assumption can be presented in the form analogous 
to Eq. (EE]): 

V|,(6) -/)w|- = 0, (27) 

where Fq is the impedance tensor of the Ith wave in the "virtual" air sphere. Hence, we 
have the relation 



r{{a) -l)^t=[P,{b) -I). (28) 

Impedance tensor F^^ contains material parameters e^^^ and /^'-^^ of the sphere inside the 
cloaking shell. At the same time, zero scattering should be held for arbitrary e^^^ and 
jj,^^^ of the inner core. It implies that the partial derivative of Eq. fl28|) with respect to 
e^^^ needs to be zero for arbitrary e^^^ to satisfy the zero scattering condition. Note that 
only the impedance tensor F'^^ contains e^^\ and therefore the right-hand side of Eq. fl25]) 

vanishes after the differentiation, which results in — ^ — —Q^; + ( r[{a) —I ) = 

l^i- ) fife = -^-(T) [ I ) fifo = 0. It is now straightforward that we need the relation to 



be satisfied, i.e., ( / ) = 0. 



By multiplying this equation by T[{a) and subtracting it from Eq. fl28|) . we arrive at the 
equation: ^0 / ^ = ^ — Fq(6) / ^, which does not contain the material parameters of 
the inner sphere. Finally, we derive the evolution operator of the cloaking layer: 

^t=i : , " I • (29) 
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This condition defines the cloak and can be satisfied for the specially chosen evolution 
operator of the cloaking shell. The evolution operator is the degenerate block matrix, 
whose inverse matrix is not defined. It should be noted that relation Eq. fl29l) is indepen- 
dent on the material of the inner sphere. On the other hand, the derived relation connects 
the wave solutions in the cloak (evolution operator Q^) with wave solutions in the "equiva- 
lent" homogeneous air sphere (impedance tensor Fq). Therefore, it effectively performs the 
coordinate transformation for the solutions, but not for the material parameters as usual. 

Substitution of the cloaking condition H'™(6) = into Eq. (12T]) results in 



From Eq. ([29D, it becomes clear that li^{^{a) = and E^^{a) = r{{a)H^{^{a) = 0. Thus one 
may conclude that both electric and magnetic fields equal zero at the boundary r = a, and 
therefore there is no field at any spatial point inside the inner sphere. In the cloaking shell, 
the fields equal zero at the inner boundary r = a (owing to the continuity of the tangential 
fields) and equal incident fields at the outer boundary r = b. Then, the field inside the cloak 
must obey (see Eq. (!20|) ) 



III. NONLINEAR TRANSFORMATION BASED SPHERICAL CLOAKS 

Now, let us consider a novel class of the nonlinear transformation based (NTB) spherical 
cloak, whose EM interaction can be characterized by the proposed scattering theorem. Fig. [1] 
can be regarded as the compressed space (r) from the original space (r'), i.e., an air sphere 
< r' < 6. We propose a nonlinear transformation function 



which obviously satisfies the transformation (when r = a, r' = 0; and when r = b, r' = b). 
The value of "x" is a factor to control the nonlinearity degree in the transformation, which 
can be arbitrary from to oo. 

Due to the invariance of Maxwell's equation under coordinate transformations from the 
original space to transformed space, the parameters (e, Jl) in the shell of Fig. [1] can be 




(30) 




(31) 




(32) 
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expressed in terms of those parameters in the original space, i.e., s' = 1 and ^t' = 1, 

£ = AA^I det ( A) , 71 = AA^I det {A) , (33) 

where A is the Jacobian matrix with elements Aij = dri/dr'-. 

One can see that the proposed prescribed function Eq. fl32|) is only dependent on radial 
position r. Then it is easy to find that the Jacobian matrix is diagonal, and Eq. (|33l) can 
thus be rewritten as 

£ = 71 = diag[A^, A^, Ajl/Ar-AeA^ = diag[-^, -^], (34) 

M'^cj) '^r^(f> \M 

where the principal stretches of the Jacobian matrix are 

^ dr {b — aYr^^^ ^ ^ r [b — aYr^^^ 
^ dr' xa6^+^(r — a)^^^ ' ^ j-i i)X+i(j, _ 

Finally, one can obtain the parameters of the NTB cloak (a < r < 6) in Fig. [T] 

6^+^(r - a)^+^ 



fir 



xa{b — a)^r 



eg= fi9 = e^ = _ ^y^^^+i ■ (36) 

Such ideal NTB spherical cloak is difficult to be fabricated in practice. However, to some 
extent, it can be alleviated by dividing the inhomogeneous cloak shell into homogeneous 
multilayers. The case for cylindrical cloaks has been studied and it shows that only several 



optimized layers can achieve the invisibility 



331]. Here, the optimization is out of the scope 



of this paper. Our paper is to reveal some novel NTB spherical cloaks which provides 
better invisibility performance than Pendry's classic one, based on the proposed general 
scattering theory. The realistic NTB cloaks can be produced using sputtering techniques, so 
that a number of discrete layers should be applied over the spherical core. To demonstrate 
the capability of the proposed spherical cloaks, we present the differential cross-sections 
normalized by the geometric cross-section of the cloak (see Eq. (1251) ). 

In Fig. [21 we analyze the dependence on the total number of the layers dividing the 
cloaking shell. The increase in the layer number gives rise to more accurate approximation of 
the original inhomogeneous model in Fig. [1], and the decrease of the scattering cross-section 
is expected with the increase of the number A^. If one uses the present scattering method 
with = 50 to divide both Pendry's linear cloak and a specific NTB cloak at x = 1, 
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the forward scattering is approximately the same, but over the whole range of scattering 
angles Pendry's cloak presents better invisibility. The following discussion will address the 
importance of this nonlinear factor x in beating the classic linear cloak. 




Scattering Angie e (degree) 

FIG. 2: Differential cross-section da /{ag sin 9 d9) of the NTB spherical cloak (x = 1) for different 
number of the layers dividing the inhomogeneous coating. Parameters: e^^^ = 1.45^, fj,^-^^ = 1.0, 
koa = vr, kQb = 2tt. 

The parameter x is a convenient tool to control the quality of the NTB spherical cloak. 
We assume the number of homogeneous sublayers = 50 for all following simulations. 
In Fig. [3] the differential cross-sections at different x are demonstrated. If x is less than 
unity, the cross-section is inversely proportional to x, which is not desired in the sense of 
invisibility. For NTB spherical cloaks with x < 1, the cloaking performance is degraded due 
to the abrupt increase of the transverse dielectric permittivity et near the inner interface 
r = a of the clad (see Fig. HI). 

The abrupt change of the material parameters is undesirable not only for invisibility 
performance but also for the realization point of view. The radial dielectric permittivity 
behaves in a similarly monotonic way for all values of x. The dependence appears to be 
mainly linear for small parameters x. The dependence of the transverse dielectric permittiv- 
ity is more complicated. One particular NTB cloak is realized at the parameter x = 2 when 
the transverse permittivity in the cloak becomes non-monotonic and eventually returns to 
Et = 2 aX r = b, which provides even lower cross-section over whole observation angles than 
the Pendry's cloak does. If we compare the dielectric permittivities of the proposed NTB 
cloak with that of Pendry's cloak (Fig. H]), it can be noted that the dependence of radial 
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Scattering Angle e (degree) 

FIG. 3: Differential cross-section of tlie nonlinear cloak with different parameters x. Parameters: 
e(i) = 1.452, = 1.0, koa = vr, kob = 2-k,N = 50. 




Clad Thickness k d 



FIG. 4: Radial Sr and transverse ej dielectric permittivities of the cloaking shell with different 
parameters x. Cloak is extended from koa = vr to kob = 2tt. 

permittivity are still quite close to each other in the cloak region. However, one may 
ask whether a; = 2 is the only choice or not. In Fig. [5], it gives the answer that in the 
sense of total cross section, there is a range of x in which the proposed NTB spherical cloak 
outperforms the classic linear spherical cloak. When "x" increases and jumps out of this 
optimal region, the cloaking effects compared with Pendry's cloak are degraded which can 
be verified by the bistatic cross-section of x = 1 and a; = 3 in Fig. [31 

Furthermore, we investigate how the transverse permittivity St varies near the optimal 
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-50-1 




-60-1 . , . , . , . ,— 

1,0 1,5 2,0 2,5 3,0 

Parameter x 

FIG. 5: Scattering cross-section versus x for NTB spherical cloaks. Parameters: e*-"^^ = 1.45^, 
/i^-*^) = 1.0, koa = vr, kob = 2tt, N = 50. The range of x where SCS is lower than that of Pendry's 
spherical cloak is the optimal region of x for desired NTB spherical cloaks. 




FIG. 6: Variance of transverse permittivity £t along the radial direction in the region of the cloaking 
shell {a < r < b) under different values of x near the optimal range as shown in Fig. [5j Parameters: 
k^a = vr and kob = 27r. 

region of "x" in Fig. [51 From Fig. [6|, it can be observed that: 1) when x is slightly above 1, 
the requirement of 6t near the inner boundary r = a drops significantly compared to those 
curves whose "x" is smaller than 1 in Fig. \5\ 2) when "x" becomes larger and larger, the et 
at the outer boundary r = b turns to be more deviated from that value of its corresponding 
Er', 3) when x falls into the optimal region, those curves of transverse permittivities are 
non-monotonic along the radial direction in the cloaking shell, and their maxima and overall 
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values of et are smaller than those whose x becomes further smaller or larger. These explain 
why there exist an optimal region for "x" where the total scattering cross-section can be 
lower than Pendry's classic one. Also, it provides us another way to predict whether a 
certain "x" for a NTB spherical cloak is optimal or not. 

Now we continue to study the dependence of its invisibility upon the ratio b/a of the 
particular NTB spherical cloak with x = 2 which is discretized into = 50 layers. We keep 
the inner radius a unchanged. In Fig. El^a), different ratios of b/a are considered. Compared 
with the other three values of b/a, it seems that b/a = 2 provides the best cloaking effects at 
nearly all angles (except for the angle at 52°) for the x = 2 NTB cloak. Another interesting 
finding is that: when b/a > 2, the cross-section will be higher than that of 6/a = 2 over the 
whole range of angles; when 6/a — > 1, though the angle-averaged cross-section will still be 
higher but at certain angles, its cross-section could be lower than that of b/a = 2. From the 
view of total scattering, it is important to consider how the ratio b/a should be selected for 
X = 2 NTB cloak so as to provide improved cloaking. In Fig. [7](b), one can clearly see the 
optimal domain of b/a in which the cross-section is smaller than Pendry's spherical cloak. 
Certainly, all values x within the desired region for the purpose of improved cloaking in 
Fig. [5] have their corresponding domain of optimal b/a. 




20 40 60 80 100 120 140 160 180 1,5 2,0 2,5 

Scattering Angle e (degree) p^gd^ ^/^ 



FIG. 7: The role of core-shell ratio b/a in the cloaking improvement for x = 2 NTB spherical cloak: 
(a) differential cross-section versus angle at selected ratios; (b) scattering cross-section versus ratio 
b/a. Parameters: e^^^ = 1.45^, ^^^^ = 1.0, koa = vr, and N = 50. 

Let us compare the scattering diagrams for both a bare glass sphere and a cloaked glass 
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sphere. We take the clad in the form of the nonhnear cloak with parameter x = 2. From 
Fig. [HI we see that the cloaking shell noticeably reduces the scattering. Another general 
property of the cloaks, i.e, the exactly diminished backscattering, is also present in the 
figure. In what follows, we consider their respective near-field wave interactions, which 
correspond to far-field results in Fig. [HI 



FIG. 8: Differential cross-section of a glass sphere with and without cloaking shell. Parameters: 
e(i) = 1.452, ^(1) = 1.0, koa = vr, kob = 2tt, x = 2, N = 50. 

The near-field perturbation of the cloaked and non-cloaked glass particles is demonstrated 
in Fig.[9l Comparing Fig.[9]^a) with Fig.[9](b), the invisibility performance is well pronounced. 
In Fig. [9](a) the EM wave travels only through the clad and takes near zero values in the 
vicinity of the inner radius a. The field does not enter the glass core, therefore it does not 
matter that which material is situated therein. If the cloak is less ideal than that shown in 
the figure, the incident field will be scattered by the spherical particle and will penetrate 
the glass core, i.e., the object becomes visible. 

IV. CONCLUSION 

We have studied the manifold of the nonlinear cloaks differing in parameter x. Since 
there is no closed form solution for the proposed nonlinear cloaks, an approximate model by 
replacing an inhomogeneous shell with homogeneous spherical layers has been numerically 
analyzed, with the help of the proposed scattering algorithm for multilayered rotationally 
anisotropic shells. The general cloaking condition was derived from the scattering algorithm, 
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FIG. 9: Real part of the electric field on x-z plane scattered (a) by the cloaking shell gathered 
round the glass core and (b) by the glass core itself. Parameters: e^^^ = 1.45^, n^^^ = 1.0, koa = tt, 
kob = 27r, X = 2, iV = 50. 

which is in contrast to the method of coordinate transform. We have also demonstrated that 
better approximate spherical cloaks can be realized by properly choosing parameters x and 
b/a. In practical applications, such a class of NTB spherical cloaks can provide improved 
invisibility performance. 

Appendix A 



In order to exclude the constant vector Hi , one should multiply Eq. (EI]) by fi^ = {^1) 



b\-l. 



r'(6) ) I T[ia) 



Hf (a) 



(A-1) 



Then Eq. flA-ll) is further multiplied by the block matrix | T[{a) —I ), and the right-hand 
side vanishes: 

/ 



T[{a) -/)^^^W|-+(rUa) 



iiim = 0. 



(A-2) 



On the other hand, vector H^™ can be obtained in a similar way. 
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